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Abstract 

We compute the energy spectra of the gravitational signals emitted when a 

mass rriQ is scattered by the gravitational field of a star of mass M » tuq. 

We show that, unlike black holes in similar processes, the quasi-normal modes 

of the star are excited, and that the amount of energy emitted in these modes 

depends on how close the exciting mass can get to the star. 
PACS numbers: 04.30.-w, 04.40.Dg 
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I. INTRODUCTION 



In this paper we study the gravitational signals emitted when a mass m is scattered by 
the gravitational field of a star of mass M. Under the assumption that itlq is much smaller 
than M, this study can be done in the framework of a first order perturbation theory, 
by assuming that the perturbations of the gravitational field and of the fluid composing the 
star are excited by the stress-energy tensor of a point-like scattered mass (a particle). 

Since 1971, when M.Davis, R.Ruffini, W.H. Press and R.H.Price JTJ computed the energy 
and the waveform emitted when a particle is radially captured by a Schwarzschild black 
hole, the perturbations of rotating and nonrotating black holes excited by infalling, scat- 
tered, or orbiting masses, have been extensively investigated. No comparable attention has 
been dedicated to the study of stellar perturbations excited by small masses, although this 
problem is certainly interesting. Indeed, we know that neutron stars exist, and observa- 
tions allow to infer the location and the mass of many of them, as well as their rotation 
periods. In addition, the worldwide experimental effort spent to detect gravitational waves, 
will hopefully be crowned by success in a not too far future. Therefore, it is time to gain 
new theoretical insight on the characteristics of the gravitational signals emitted by neutron 
stars. Until now, the application of the theory of stellar perturbations has mainly regarded 
the determination of the frequencies of the quasi-normal modes, at which stars are expected 
to pulsate and emit gravitational waves. However, whether these modes can be excited, and 
how much energy is emitted at the corresponding frequencies, can be understood only by 
considering realistic astrophysical situations. The theory of stellar perturbations, suitably 
matched with the theory of black hole perturbations, allows to explore this new field. 

In our investigation we shall consider a polytropic star with mass and radius typical of 
neutron stars, and integrate the equations describing the scattering of a small mass. In the 
interior of the star, we shall use the equations for the axial and polar perturbations derived by 
S.Chandrasekhar & V.Ferrari in 1990 in ref. ||, (to be referred to hereafter as Paper I), that 
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hold also in a more general gauge appropriate to describe non-axisymmetric perturbations. 
Outside the star, one may continue the solution by integrating the Regge- Wheeler and the 
Zerilli equations 0,0, with a source term given by the stress energy tensor of the scattered 
mass moving along a geodesic of the unperturbed spacetime. However, it is known that the 
source term of these equations diverges when the mass reaches the periastron, and therefore 
we switch to a different formalism, and use the generalized non-homogeneous Regge- Wheeler 
equation, which was introduced by T.Nakamura and M.Sasaki to overcome this problem. 
We shall integrate this equation by adopting the procedure discussed in ref. f| (to be referred 
to hereafter as Paper II), and further developed in this paper. We shall show how the problem 
of matching the interior and exterior solution can be solved, and find the energy spectra of 
the emitted radiation for different values of the orbital parameters of the scattered mass. 

The plan of the paper is the following. In section II we shall briefly review the equations 
to be integrated in the interior of the star; in section III the equations governing the exterior 
perturbations will be shown; the procedure to find the complete solution and the matching 
conditions will be discussed in section IV; in section V we will describe a model of polytropic 
star for which we shall compute the characteristic oscillation frequencies; in section VI we 
will show the energy spectra emitted when a small mass is scattered by such star. 

II. THE EQUATIONS DESCRIBING THE PERTURBED SPACETIME INSIDE 

THE STAR. 

In order to describe the perturbations induced by a mass scattered by a spherical, non- 
rotating star, we choose a gauge appropriate to describe non-axisymmetric perturbations 
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and Yimid, p) are the scalar spherical harmonics. The perturbed part of the metric ( ^.1| ) 
has been Fourier-expanded, and the usual decomposition in tensor spherical harmonics has 



been performed (cfr. HU). The functions N im (r,u), L im (r,uj), 



are 



the radial part of the polar, {even) metric components, whereas h® m (r, lu) and hj m (r,u) 
are the axial (odd) part. It should be mentioned that in the axisymmetric case (m = 0) 
the previous gauge reduces to that used in Paper I for the polar perturbations, and to the 
Regge- Wheeler gauge || for the axial ones. 

The unperturbed metric functions u(r) and /i2( r ) have to be determined by numeri- 
cally integrating the Einstein equations coupled to the equations of hydrostatic equilibrium, 
for an assigned equation of state. Under the assumption that the star is composed by a 



perfect fluid with energy-momentum tensor given by T Q/3 = (p + e] 



u u 
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and e are, respectively, the pressure and the energy density, the relevant equations are 
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where m(r) = J r er 2 dr is the mass contained within a sphere of radius r. The solution for 
the function v(r) requires the determination of the constant z/ , which can be found by 
imposing that at the boundary of the star the metric reduces to the Schwarzschild metric 



-_ R = 1- 2M/R, 



(2.4) 



where M = m(R) is the total mass. 

It is easy to check that, since Y im (fl, tp) = (— \) m Yl_ m ($, ip), the perturbed metric 
functions satisfy the following property 
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where F em (r,uj) is any of the functions N em (r, u), L im (r, u), V em (r, u), T em (r, u) , and 

By writing explicitely the perturbed Einstein's equations coupled to the hydrodynamical 
equations in the interior of the star, it is possible to verify that the resulting separated 
equations coincide with those given in Paper I, both for the polar and for the axial per- 
turbations. As a consequence, the decoupling of the gravitational perturbations from the 
perturbations of the thermodynamical variables that was performed in Paper I for the polar 
equations, is possible also in the present case of non-axisymmetric perturbations. Therefore, 
the equations to integrate inside the star, to find the values of the polar and axial functions 
at the boundary, are the following. If we consider stars with a barotropic equation of state, 
the polar equations are (cfr. eqs. (72)- (75) of Paper I) 
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where Q = e, r /p,r, the function V has been replaced by X = nV, with n — (£— 1)(£+2)/2, 
and the function G takes the place of T 
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Outside the star, where e and p vanish, these equations become appropriate to describe 
the Schwarzschild perturbations, and can be reduced to the Zerilli equation [|J 
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The value of the Zerilli function at the surface, Z pol (uj,R), can be found in terms of the 



solution of the eqs. ( |2.6|) as follows (cfr. eq. (93) Paper I) 
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and similarly for its first derivative. For simplicity, in eqs. ( |2.6| - |2.9|) we have omitted the 
harmonic indices (£m). 

The equations for the axial perturbations can be reduced to a single wave equation (cfr. 
eqs. (148)and (149), Paper I) 
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where r* = J r e u+ ^ 2 dr, and the function is related to the axial metric components 

by 
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Outside the star eq. (|2.10|) reduces to the Regge- Wheeler equation. 



III. THE EQUATIONS DESCRIBING THE PERTURBED SPACETIME 

OUTSIDE THE STAR. 

The metric exterior to a nonrotating star reduces to the Schwarzschild metric, and there- 
fore the perturbed spacetime is described by the perturbations of a Schwarzschild black 
hole. When the source exciting the perturbations is a scattered mass, the source term of 
the Regge- Wheeler and of the Zerilli equations diverge at the periastron (where j- = 0), 
and alternative equations have to be used. In 1973 J.M. Bardeen and W.H. Press derived 
a wave equation with a complex potential barrier for the Weyl scalars and r 4 5^4 
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0. This equation was subsequently generalized by S.A. Teukolsky ||, who found a master 
equation governing the electromagnetic, gravitational and neutrino-field perturbations of a 
Kerr black hole, which reduces to the Bardeen-Press equation when the rotation parameter 
is set to zero. We shall refer to the Bardeen-Press- Teukolsky equation as the BPT equation. 
After separating the variables, the radial BPT equation can be written as 

"(rV 2 + Ai(r - M)r 2 u) 
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where A = r 2 — 2Mr, the BPT function is related to by 
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and _ 2 5 , £ m ( , (9, if), in terms of which the separation is accomplished, is the spin-weighted 
spherical harmonic 
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The source Ti m (uj,r) is given by 
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and the factors T^ M (r) are the radial parts of Newman- Penrose components of the stress- 
energy tensor (see Paper II, app. A. 2). In the present paper we prefer to use another 
equation, related to the BPT equation, which was derived by Nakamura and Sasaki ||. 
They showed that, given a solution of the inhomogeneous BPT equation , there exists a 
function Zj^(u,r) related to ^f^iu, r) by the Chandrasekhar transformation M 
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where 
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which satisfies the generalized inhomogeneous Regge- Wheeler equation 
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The source term is related to the source of the BPT equation by 
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We shall now write the source term of the Nakamur a- Sasaki equation explicitely. If we use, 
as in paper I, the convenction = 2T fll/ , the stress-energy tensor of a point-like mass 

mo, moving along a geodesic r(r) of the unperturbed spacetime, is 
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where Tj(r) are the solution of r(r) = f, the solution can be divided in two branches, 
corresponding to the incoming and outgoing part of the trajectory. Consequently, the stress- 
energy tensor can be written as 

1 2 dx* dx v 
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where U(r), Qi(r) are the time and angular position of the particle on the i— eth branch of 
the trajectory. In terms of the tensor ( [3.11| ), the source of the BPT equation can be written 
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the prime indicates differentiation with respect to r, and _o^ m (^( r )) >-i ^ m (^( r )) are 
the complex conjugate of the spin- weighted sperical harmonics of weight and 1, and a = 



sign 



dr 
dr 



. We shall assume that a particle mo starts its journey at radial infinity with 
energy E and angular momentum L z , and we shall put <^(r t ) = r(r t ) = 0, at the turning 
point r = r t . We choose E and L z , such that the particle is scattered by the star, and 
follows a geodesic on the plane ■& = ~ , described by the equations 
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The source term of the Nakamur a- Sasaki equation can be derived from eqs. (|3 . 9| ) , ( pTT2|) and 
fl3.13j ), as in |15|]. The result is 
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The constants K ,... are 
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IV. THE SOLUTION OF THE NAKAMURA-SASAKI EQUATION AND THE 

MATCHING CONDITIONS. 



In order to integrate the Nakamura-Sasaki equation (|3.8|) outside the star, we need to 



know the value of ZfJ^ and its derivative at the surface r = R. The Nakamura-Sasaki 
function is related to the solution of the Regge- Wheeler equation for the axial perturbations, 
Z ax , and to that of the Zerilli equation for the polar ones, Z po1 , by the following relations 
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It should be stressed that both Zj m (u,r) and Zf m (uj,r) satisfy the Regge- Wheeler 
equation. Thus, Z^(u,R) and its derivative can be found in terms of Z pol (u,R) and 
Z ax (u, R) and their first derivatives. These values can be found by numerically integrating 
the equations for the polar and axial perturbations in the interior of the star (cfr. eqs. [T9] and 
2.10| ), by imposing the boundary conditions that, for each assigned value of the frequency, 
the solution is regular near the origin, and that the perturbation of the thermodynamical 
variables vanish at the boundary, so that the perturbed spacetime reduces to vacuum (cfr. 
Paper I, sec. 6,7 and 11). However, these conditions do not define the amplitude of the 
functions Z pol (u),R) and Z ax (uj,R), which depends on the exciting source. If we name 
Z^(uj, R) and Z^(u, R) the values found by integrating the interior equations with an 
arbitrary amplitude at the centre of the star, the 'true' values of Z pol (u,R) and Z ax (uj,R) 
will be 
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and xlmi. 00 ) an d Ximi^) have to be found. Consequently, the functions Z} m (uj,R) and 
Zj m {uj,R) defined in eqs. ( |4.2| ) suffer the same ambiguity 
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Zemi^iR) = xl m (u) Z$ m (u,R) 
and Z\ m {uj , R) , Zf m (uj , R) are (cfr. eqs. [4] 
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x\ m {oj) and x\ m {uS) determine how much of the polar and of the axial modes are excited 
by the particular source we are considering. 
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Eq. ( p.8|) has to be solved by imposing that, at the surface, the solution of the perturbed 
equations in the interior matches continuously with the exterior solution, which, in addition, 
has to behave as a pure outgoing wave at radial infinity, i.e. 
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The solution of eqs. ( |4.8|) can be found by defining the Green's function appropriate to the 
problem, as explained in Paper II. Here we briefly summarize the procedure. For r > R 
we define the function Y^lu, r*), which satisfies the homogeneous Regge- Wheeler equation 

with the pure outgoing wave condition at infinity 

L RW Y e (u,r*) = 

(4.10) 

Y t (u, oo) = e iujr * 

and the functions Z\{uj,r^) and Z|(u;,r*) that, as Y^u, r*), satisfy a homogeneous 
Regge- Wheeler equation, with boundary condition at the surface of the star given by eqs. 
( [4.6| ) and ( f4.7|) , respectively. It should be noted that, since the Regge- Wheeler operator does 
not depend on m , the functions Yi(uj,r*), Zj(uj,r*) and Z]{uj,r^) will be independent 
of m as well. We then construct the solution of the Nakamura-Sasaki equation as follows 
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where Wi(uj) and WfiuS) are the wronskians 
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and a£ m {uj) are constants. It should be mentioned that, since the source S^(uj, r*) 



given in eqs. (|3.14|) -( [3~T7|) diverges as 7 1 at the turning point, the actual integration of 



eq. ( 4.11 ) can be performed by switching to the proper time r(r). 



It is easy to check that the solution (|4.11| ) satisfies the pure outgoing wave condition at 
infinity, and that the matching condition at the boundary are fulfilled provided the constants 
xjmi. 00 ) an d Xemi^) are defined as 
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The solution at radial infinity, to which we are primarily interested, therefore is 
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The constants ac£ m (uj) will be determined by the following arguments. Due to the properties 
of the spherical harmonics 

Y tm (#, ip) = (-l) m Y em ($, ip + 7t) = {-l)%m{K - V + tt) = (-l) £+m F, m (vr -$,<p), (4.16) 

by an inspection of the behaviour of the source S^- , it is easy to verify that the Nakamura- 
Sasaki function satisfies the following identity 
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and consequently 
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In addition, by explicitely evaluating 5^/^(t, r, ip) for the perturbed metric ( |2.1|) , and by 
using eq. (|3.6|), we find 
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From the property of the metric perturbations given in eq. ( |2.5| ), and since D* ch (uj) = 
Dch(-u) it follows that 

Z^(u,r) = (-l) m Zl m (-u,r) (4.20) 
^( W| r) = (-ir^Ur). 
Eqs. ( |4.18[ ) and (4.2U) are compatible only if 

{-\) l Z\ m {u, R) = (-l) m ZU", R) (4.21) 
-(-lYZ] m (u;,R) = (-irZ* m (u,R), 

which imply that 

• if (£ + m) is even, Z] m (u,R) = 0, 

• if (£ + m) is odd, Z} m (u,R) = Q, 

From eqs. ( (4.5|) we see that if (£ + m) is even, in order Zf m (u, R) to be zero x 2 m (u;) 
must vanish, which means that ac£ m (uj) = 1 (see eq. |4.14p . Similarly, if (£ + m) is odd 
Xe m ( UJ ) must vanish, i.e. ae m (u>) = 0. Thus, the complete solution of the Nakamur a- Sasaki 
equation at radial infinity is (see eq. |4.15|) 
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Equations (|4.22 ) and (|4.23| ) show that the polar and the axial perturbations contribute to 
the solution Zf^ depending on {I + m) being even or odd, respectively. This fact is 
not surprising. For example, if we consider the non-homogeneous Regge- Wheeler equation 
written for the process under consideration 
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(4.25) 



L/), n(£) = W+l)}~\ and m(f) = [£(£+!) (£+2) (£-l))-\ 



when (£+m) is even Yl m ^ \$ =7T / 2 = 0, the source term is zero, and the axial perturbations 
are not excited. 

It should be mentioned that the procedure used in this paper to find the constants 
ctimiu) exploits the symmetry properties of the source of the Nakamura-Sasaki equation, 
and is much simpler than that described in Paper II. 
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V. THE MODEL OF STAR 



We have integrated eqs. ( |2.6|) and ( p. 10 ) in the interior of a star with a polytropic 
equation of state p = Kp n , with K = 100 km and n = 2. If the central density 
is chosen to be p c = 3 • 10 15 g/cm 3 , the radius and mass of the star are, respectively, 
R = 8.86 km, and M = 1.266 M , with a ratio R/M = 4.7 For this model of 
star, in Fig. 1 we plot the resonance curve evaluated for the £ = 2, polar perturbations, 
which allows to locate the frequencies of the quasi-normal modes of the star. This curve 
is obtained by integrating the polar equations inside and outside the star (eqs. [2T^ , p.9[ ) for 
assigned real values of the frequency, and by evaluating the amplitude of the standing wave, 
[a 2 (o;) + (3 2 {uo)} , prevailing at radial infinity. It can be shown that the values of frequency 
at which this curve exhibits a sharp minimum correspond to the real part of the complex 
eigenfrequency of a quasi-normal mode, provided the corresponding imaginary part is small 
enough {oji « u> .) @, (T^ . The damping time associated to a mode is related to the 
curvature of the parabola that fits the curve near a minimum: smoother minima correspond 
to shorter damping times. 

As in newtonian theory, the classification of the polar modes is based on the behaviour 
of the perturbed fluid according to the restoring force that is prevailing []n}| : the g- modes, 
or gravity modes, when the force is due to the eulerian change in the density, the p-modes, 
when it is due to a change in pressure, and the f-mode, that is the generalization of the only 



possible mode of oscillation of an incompressible homogeneous sphere ||12|| . An inspection 
of the behaviour of the thermodynamical variables in correspondence of the minima of the 
resonance curve shown in Fig. 1, allows to identify the corresponding modes. The resonance 
curves for i = 3 and I = 4 are plotted in Fig. 2. A resonance curve can be computed 
also for the axial perturbations by solving eq. ( |2.10j ). For the model of star we consider, it 
exhibits a monotonically decreasing behaviour, showing that no slowly damped axial modes 
exist in this case. The algorithm based on the resonance curve allows to find only the slowly 
damped modes « u .). To find the highly-damped modes (the w- modes (T^|) other 
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methods have to be used. In table 1, we tabulate the values of the complex frequency of the 



first few polar w-modes computed in ref. |]T3[ for the same model of star considered in the 
present paper. We do not go into more details about these modes because, as we shall see, 
they are not excited in scattering processes. 



VI. NUMERICAL RESULTS 



As discussed in section IV, we find Z^ o1 (uj,R) and Z^ x (uj,R) and their first derivatives 
by solving the interior equations ( |2.6|) and ( |2.10|) for assigned values of the frequency. With 
these values we construct Z}(u, R) and Zf(u, R), and integrate the homogeneous Regge- 
Wheeler equation for these functions and for the function Yg(u},r*) defined in eq. ( [4.10|) , 
up to radial infinity. These functions are needed to construct the solution of the Nakamura- 
Sasaki equation, and find its asymptotic behaviour given by eqs. ( |4 . 2 2| ) and ( |4.23| ) . The 
perturbations are assumed to be excited by a massive particle mo, which is scattered by 
the gravitational field of the star; we set the angular momentum and the energy respectively 
to L z = 5, and E = 1.007. For these values, the turning point is located at r t /M = 9.2 
(case a). We then repeat the calculations allowing the particle to get closer to the star: we 
choose E = 1.097, so that r t /M = 5.0 (case b). 

For a plane wave (see for example |L4j] pg. 522 eq. (488)) 



1 r-. 



h w (t, r, 0, 6) + ih^p(t, r, 0, 6) 



(6.1) 



and since the (Or)-component of the stress-energy pseudotensor can be written as 
( dE\ 1 
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(6.2) 



it follows that the energy emitted in gravitational waves per unit solid angle and unit fre- 
quency is 

dE 



dVLdoj 2u 2 



|5v|> 4 ( w ,r,^)| 2 . 



(6.3) 
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By integrating over the solid angle, and by using the relation existing between 6^4 and 
the Nakamur a- Sasaki function and discussed in section III, it can be shown that the energy 
emitted per unit frequency can be expressed as a function of the amplitude of Z$~ at 
infinity, i.e. 

dE lm _ 2 NS 2 



16u/|0)| 2 , 0<cu<oc, (6.4) 



where 



^>,°°) - C(^)e w *, (6.5) 

and Af^iyj) is determined by numerical integration. We find that, as in the case of 
the scattering of masses by a black hole, the energy spectrum of the emitted gravitational 
radiation is contributed mainly by the i = m component. 

In Fig. 3 we plot the I = m = 2 energy spectrum emitted in case a, when the turning 
point is r t = 9.2 M. It is interesting to compare this spectrum with that obtained when a 
black hole scatters a small mass (see \TE\ for an extensive review). In that case the energy is 



emitted essentially by the scattered mass, and most of it is radiated when the mass transits 
through the turning point. Indeed, the energy spectrum is peaked at a frequency which is 
related to the angular velocity of the mass at the turning point 

<*. = l(%) ■ (6-6) 



v dt > 



r=r t 



Thus, the black hole quasi-normal modes are not excited in these processes. In our case a, 
and for I = 2, the frequency corresponding to the angular velocity of the mass at the 
periastron is u n M = 0.092 and the spectrum shown in Fig. 3 exhibits a sort of peak 
close to that frequency, showing that part of the energy is still emitted by the mass as a 
synchrotron radiation. However the very distinctive feature of that spectrum is a very sharp 
peak which occurs at the frequency of the fundamental f-mode (cfr. Fig. 1). If we allow 
the mass to get closer to the star, as shown in Fig. 4 for case b when the turning point is 
r t = 5, the amplitude of the f-mode peak increases by more than a factor 10 2 , and new 



peaks appear, corresponding to the excitation of the p-modes. (In this case u> n M = 0.221). 
Thus, we can conclude that, unlike the case of black holes, the polar quasi-normal modes of 
the star can be excited in scattering processes, to an extent that depends on how close the 
mass can get to the star. Conversely, from table 1 and Fig. 3 and 4, we see that the w-modes 
are not significantly excited in these processes. It should be mentioned that, although we 
plot the energy spectra for uoM < 1, we have extended our numerical integration at higher 
frequencies finding that the signal is negligible compared with that shown in the figures and 
no further modes appear. In Fig. 5 we plot the energy spectrum for £ = 2 and different 
values of — £ < m < £. We see that when £ + m is odd, no modes are excited, because 
the spectrum is contributed only by the axial perturbations which are not resonant for this 
model of star. In Fig. 6 we compare the energy spectrum emitted for £ = 2,3 and 4. The 
excitation of the quasi-normal modes is exhibited by all multipoles. 



VII. CONCLUDING REMARKS 

In 1983 Lindblom and Detweiler investigated the relation existing between the frequencies 
of oscillation of stars and the equation of state (EOS) prevailing in the interior, showing 
that quasi- normal modes carry information on the internal structure of a star jig] . They 
computed the frequency of the fundamental mode of oscillation (the f-mode), for several 
EOS's, and more recently this study has been extended to include a few modern EOS's 
suggested for neutron stars, for which also the frequencies of the first p-mode and of the 
polar and axial w-modes have been computed. Besides, it has been shown that the knowledge 
of these frequencies allows to infer empirical relations between the mode frequency and the 
macroscopic parameters of the star: the mass and the radius [fLTfl , [ |i~8|j . Thus, there is no 
doubt that the observation of the quasi-normal modes would provide relevant information 
on both the structure of neutron stars and on the nature of the nuclear interactions at 
supernuclear densities. 

The real problem is to ascertain whether these modes can be excited in realistic astro- 
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physical situations, and what is the amount of gravitational energy that is emitted at the 
corresponding frequencies. The work presented in this paper is a first step in this direc- 
tion. The energy spectra plotted in Fig. 3-6 clearly show that the quasi-normal modes of 
stars can be excited in scattering processes. The waveforms of the corresponding signals 
are essentially exponentially damped, pure sinusoids, at the frequency of the f-mode of the 
considered star ( v ~ 3 kHz,r ~ 0.07 s), and with an amplitude which scales as itiq 2 , 
and depends on how close the scattered mass gets to the star. 

The p- modes can be excited if the mass gets sufficiently close (Fig. 4), but the f-mode 
contribution is likely to be the dominant one. Thus, the signal emitted by a star excited by 
a scattered mass is a pure note, emitted at the frequency, and with the damping time, of the 
fundamental f-mode. The spectra we compute do not show the excitation of the w-modes. 
However, they may be excited in other processess, like for example the capture of masses by 
stars. 

We plan to extend the present work, to compute the energy and the waveforms emitted 
by masses orbiting around compact stars. 
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FIGURES 

FIG. 1. The £ = 2- polar resonance curve is plotted versus the normalized frequency wM. 
Each minimum corresponds to a quasi-normal mode, and the modes are labelled according to the 
Cowling classification. 

FIG. 2. The £ = 3 (right) and £ = 4 (left) polar resonance curves are plotted as a function 
of the frequency as in fig. 1. 

FIG. 3. The i = m = 2 - component of the energy spectrum of the gravitational radiation 
emitted when a mass tjiq, with L z = 5 and E = 1.007, is scattered by a polytropic star, is 
plotted versus the normalized frequency. In this case the turning point is located at r t /M = 9.2 
(case a). 

FIG. 4. The i = m = 2 - component of the energy spectrum emitted when the scattered 
mass mo has an energy at infinity, E = 1.097, higher than that considered in fig. 1. In this 
case the turning point is located at r t /M = 5.0, and the mass gets closer to the star (case b). 

FIG. 5. The m-components of the energy spectrum are plotted for I = 2 for case b. We see 
that when (£ + m) is even, the polar quasi-normal modes are excited. When (£ + m) is odd, 
the energy is emitted essentially by the scattered mass (see text). 

FIG. 6. The energy spectrum emitted at £ = 2 is compared with that emitted by higher 
multipoles (case b). 
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TABLES 

TABLE I. The values of the real and imaginary part of the frequencies of the first few polar 
w-modes for the model of star under consideration. 



N Re(toM) Im(uM) 

1 0.53 0.30 

2 0.88 0.38 

3 1.22 0.43 

4 1.57 0.47 

5 1.91 0.50 
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